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ABSTRACT 

In this paper we investigate the evolution of a pair of interacting planets - a Jupiter 
mass planet and a Super-Earth with the mass 5.5 M s - orbiting a Solar type star and 
embedded in a gaseous protoplanetary disc. We focus on the effects of type I and II 
orbital migrations, caused by the planet-disc interaction, leading to the Super-Earth 
capture in first order mean motion resonances by the Jupiter. The stability of the 
resulting resonant system in which the Super-Earth is on the internal orbit relatively 
to the Jupiter has been studied numerically by means of full 2D hydrodynamical 
simulations. Our main motivation is to determine the Super-Earth behaviour in the 
presence of the gas giant in the system. It has been found that the Jupiter captures the 
Super-Earth into the interior 3:2 or 4:3 mean motion resonances and the stability of 
such configurations depends on the initial planet positions and eccentricity evolution. 
If the initial separation of planet orbits is larger or close to that required for the 
exact resonance than the final outcome is the migration of the pair of planets with 
the rate similar to that of the gas giant at least for time of our simulations. Otherwise 
we observe a scattering of the Super-Earth from the disc. The evolution of planets 
immersed in the gaseous disc has been compared with their behaviour in the case of 
the classical three-body problem when the disc is absent. 

Key words: planets formation, numerical simulation, protoplanetary discs, migra- 
tion, mean motion resonances 



1 INTRODUCTION 

More than ten years have passed since the discovery of 
first extrasolar low mass planets (with masses 0.02, 4.3 and 
3.9 M ffi ) around pulsar PSR B 1257+ 12 (IWolszczan fc Fraill 
1992) an d the first hot Jupiter o rbiting 51 Pegasi, a Solar- 
type star (|Mavor fc Quelozl 19951 ). Since then, more than two 
hundred gas giants have been found around main sequence 
stars and only four low (less massive than 10 M®) mass 
planets. The first among these four, with a mass 7.5 or- 
biting the nearby star Gliese 876 has been discovered by the 
high-precision radial ve locity monitoring a t the W. M. Keck 
Observatory in Hawaii l|Rivera et al.ll2005h . The second one, 
OGLE-2005-BLG-390Lb with a mass of 5.5 M @ , has been 
found by using the gravitational microlensin g technique by 
OGL E, PLANETS and MOA collaborations (|Beaulieu et all 
2006). Finally the last two, with masses a bout 5 and 8 M® , 
has been observed around star Gliese 581 l|Udrv et al.ll2007l ) 
thanks to HARPS (High Accuracy Radial Velocity for Plan- 
etary Searcher), a very precise spectrograph located on the 
ESO telescope at La Silla in Chile. All three stars harboring 
low mass planets are M dwarfs. The planets in this mass 
range are often called Super-Earths and we will adopt this 
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name here. These few examples are just a preview of what 
is to come soon. The microlensing and the high-precision 
radial-velocity planet-search programmes can bring new dis- 
coveries at any time. Even a network of small ground-based 
telescopes is able to increase the number of known Earth-like 
planets using the Transit Timing Variation (TTV) method. 
The French-European mission COROT, which is capable to 
detect low mass planets, is already operating and collect- 
ing data. In 2009 the Kepler satellite will be launched and 
also will search for Earth-size planets. In other words, it is 
a perfect time to have a closer look at the possible planet 
configurations formed as a result of the evolution of young 
planetary systems in which different mass planets might co- 
exist with each other. Here we focus on the resonant con- 
figurations between two planets - a gas giant (we will call 
it Jupiter) and a Super-Earth, orbiting a Solar-type star. 
Our prediction of the occurrence of the commensurabilities 
in such systems may turn out to be particularly valuable in 
the case of the TTV observations, because the differences 
in the time intervals between successive transits, caused by 
planet- planet interaction s, are largest near mean-motion res- 
onance |Agol et al.ll2005l ). 

It has been recognized that resonant structures may 
form as a result of large scale orbital migration in young 
planetary systems due to tidal interactions between plan- 
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ets or protoplanets and a gaseous disc, in which the whole 
system is still embedded. The different mass objects, which 
we expect to find in forming planetary systems, will mi- 
grate with different rates. The final configurations will de- 
pend on the intricate interplay among many physical pro- 
cesses including planet-planet, disc-planet and planet-star 
interactions. One scenario is that the convergent migration 
brings the giant planets closer to each other and they can 
become locked in low order commensurability l|Brvden et al.l 
l200d ; iKlevI l2000h as it is observed at least in three mul- 
tiplanet systems (GJ 876, HD 82943 and 55 Cnc). Also, 
the low mass planets can under go the convergent migration 
and form a resonant structure ( Papaloizou fc Szuszkiewicd 
120051 ). The pulsar planets around PSR B1257+12 might be 
an outcome of such scenario. The differences in the migra- 
tion rates of low and hi gh mass planets may also lead to the 
convergent migration IjHahn fc Ward! I1996T ) . We will con- 
sider this type of convergent migration here as a mecha- 
nism for capture of a Super-Earth by a Jupiter into mean 
motion commensurabilities and we will follow the pair of 
planets further in their evolution to establish the stabil- 
ity of the i r con figuration. It has been found already by 
iThommesI (|2005h that for typical protoplanetary disc pa- 
rameters the low mass planets can be captured by mas- 
sive planet into exterior resonances. The resonant cap- 
tures play also an important role in modelling terrestrial 
planet formation and their survival in the presence of mi- 
grating gas giant. The early divergent conclusions on the 
occurrence of terre s trial planets in hot Jupiter systems 
llLevison et al.ll200ll; I Armitagd [20031 : iRavmond et alj|2005l : 
iMandell fc Sigurdsson _ 20031) have b e en clar ifi ed by most 
recent studies dFogg fc Nelson! 120051. I2007all3; IZhou et all 



l2005l : IRavmond et al-lbOOrj ; IMandell et al.ll2007l ) which pre- 
dict that terrestrial planets can grow and be retained in the 
hot- Jupiter systems both interior and exterior of the gas gi- 
ant. The most relevant feature of these investigations for the 
present study is the possibility of the existence of terrestrial 
planets on the internal orbits relative to Jupiter. Such plan- 
ets could become captured in the mean motion commensu- 
rability which is then maintained during their evolution so 
an outcome in this case will be the Jupiter and the low-mass 
planet migrating together towards the star. 

Our aim here is to consider in full detail the evolution 
of the close pair of planets: a gas giant on the external orbit 
and a Super-Earth on the internal one in the young plane- 
tary system, when both planets are embedded in a gaseous 
disc. We have performed 2D hydrodynamical simulations in 
order to determine the occurrence of the first order mean 
motion resonances in such systems and at the same time to 
examine the possible planet configurations as the outcome 
of the convergent orbital migration. 

Our paper is organized as follows. In Section 2 we sum- 
marize shortly few facts about two types of migration, which 
we will use in the present study. In Section 3 we discuss the 
Jupiter and Super-Earth system in the context of the re- 
stricted three body problem when the disc is absent. More- 
over we have performed simple N-body calculations in or- 
der to determine the Super-Earth dynamics in a presence of 
Jupiter and compare it with the classical results of celestial 
mechanics. This provides us with a well studied and under- 
stood framework in which we can present the results of our 
hydrodynamical simulations of a young planetary system, 



where planets are still embedded in a gaseous disc. In Sec- 
tion 4 we describe our hydrodynamical simulations and show 
the planetary orbit evolution together with the changes in 
the protoplanetary disc structure. The discussion and our 
conclusions are given in Section 5. 



2 SLOW AND FAST MIGRATION 

Migration due to planet-disc interaction might play an im- 
portant role in shaping up planet configurations in the plan- 
etary systems as we have already mentioned in the previous 
section. From the dynamical point of view one of the most 
important consequences of this process is an occurrence of 
planets locked in mean-motion resonances and this will be 
our main concern here. 

The migration rates for different planet masses has 
been estimated by num ber of authors, see the review by 
iPapaloizou et all l|2006l ). Their results have been illustrated 
in Fig. [T] where we plot the migration time of a planet as a 
function of its mass. There are two mass regimes which are 
of interest here, namely (0.1 - 30M ffi ) and (150 - 1500M©) 
for which in a typical protoplanetary disc we can talk about 
two different types of migration, called simply type I and 
type II respectively. 

The migration time for low mass planets embedded 
in a gaseous disc ( type I migration) has been derived by 

m 



iTanaka et all ([2002) in the form 



rj = (2.7 + 1.17) 
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Here m p is mass of the planet, r p is the distance from the 
central star M,, E p is the disc surface density, H/r and Q p 
are the disc aspect ratio and angular velocity respectively. 
The coefficient 7 depends on the disc surface density profile, 
which is expressed as E(r) oc r Assuming 7 = for the 
flat surface density distribution, E p = 2000 kg/m 2 , H/r — 
0.05 and r v — 5.2AU we have calculated the migration time 
as a function of the planet mass and draw it in Fig. [T] 

Type II migrators open a gap in the disc and their evo- 
lution is determined by the radial velocity drift in the disc 

where v is a kinematic vi scosity parameter. The migration 
time can be estimated as (|Lin fc Papaloizoulll993l ) 
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= 5.2AU and different 
,9 - 10~ 5 and 10" 6 ) ex- 



It has been shown in Fig. \T\ for 
values of v (10 -5 , 2 ■ 10~ 5 , 3 ■ 10" E 
pressed in the dimensionless units discussed in Section 14.11 
The migration time in this formulation does not depend on 
the planet mass. In drawing lines for a given viscosity pa- 
rameter we have taken into account the condition for a gap 
opening in the disc which reads 

nhL > J!l ( 4 ) 

M, r 2 p Q p V ; 

and from which it is clear that the bigger is viscosity v the 
bigger is the mass of the p lanet, which is able to open a 
gap. Recently, lEdgarl ll2007l) has argued that similarly like 
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Figure 1. The comparison between the times of migra- 
tion for planets with different masses at 5.2 AU. We have 
concentrated on classical type I and II migration cases, for 
which the theoretical and numerical evaluations agree d very 
well and their understanding is rather established dWardl 

dl2 



19971; IPapaloizou fc Larwoo J2000l; IPapaloizou fc Terquemll2006t 
Bate et alj|2003t iD'Aneelo et alj|2003h . 



in the case of type I migration also here the migration time 
depends on the planet mass and the disc surface density. 
His relation gives for our set of parameters somewhat faster 
migration for the Jupiter in comparison with the equation 
([3|, so it does not change our argument here. 

We have chosen the mass of the type I and type II mi- 
grators to be 5.5 M($ and 1 Jupiter mass respectively. Their 
locations have been marked in Fig. [T] accordingly to their 
masses and in the case of the Jupiter also to the viscos- 
ity adopted in the disc (y = 10 -5 ). It is obvious from this 
figure that we should expect a convergent migration if the 
Jovian type planet is at the external orbit, further away from 
the star, and the Super-Earth at the internal one. We have 
checked that for the reverse configuration we are getting 
the divergent migration. For the disc with the low viscos- 
ity this situation may change and the Jupiter mass planet 
would migrate slower than the Super-Earth. In order to get 
the convergent migration the Super-Earth should than be 
lo cated at th e exter ior orbit exactly as the Earth considered 
m iThommesI (|2005h . The overall picture is quite simple and 
straightforward. For a given disc profile and masses of the 
planets, the relative speed of type I and type II migrators is 
controlled by the viscosity in the disc. However, the details 
of this scenario require more attention. We will come back 
to this issue in Section [S] 

Starting from the configuration illustrated in Fig. [T]we 
have investigated the possible resonances in the system con- 
taining a star and two planets with masses, which are much 
different from each other. Their mass ratio to the central 
star is equal to msE = 1-65 • 10" 5 for the Super-Earth and 
mj — 1 ■ 10 -3 for the Jupiter, so their masses differ from 
each other by almost two orders of magnitude and the Super- 
Earth, in the absence of the disc, could be treated at least 



approximately in the context of the classical restricted three- 
body problem, where one body is so light that it does not 
influence the motion of the two others. In the next section 
we compare the dynamics of our system with the classical 
results of the three-body problem. 



3 CLASSICAL THREE-BODY PROBLEM 

Before tackling the full problem of the Jupiter and Super- 
Earth, orbiting around Solar-type star, embedded in a 
gaseous disc let us first consider a simpler well studied sit- 
uation in celestial mechanics when the disc is not present 
in the system. We are interested in motion of three celes- 
tial bodies under their mutual gravitational attraction. In 
our system the mass of two planets is much smaller than 
the central star mass and the mass of the Super-Earth is 
significantly smaller than that of the Jupiter. 

The dynamics of a pair of interactin g planets orbitin g 
a star has been studied in full details by iGladmanl (|l993f) . 
The main conclusion of his paper is that for initially circular 
planetary orbits with the semi-major axes a and o(l+A), the 
relative orbital separation A between the two planets with 
the mass ratios to the mass of the central star equal mi and 
mi, respectively, has to satisfy the following criterion 



A > 2.4(mi + m 2 ) 



1/3 



(•») 



in order to be Hill stable. This means that for the plane- 
tary orbit separations fulfilling this condition planets will 
not suffer from the close encounters at any time. The close 
encounter will result in the scattering of the Super-Earth, 
which may leave the system, fall onto the star or the Jupiter. 
In our case, for the Jupiter (mi = mj) and Super-Earth 
(m,2 = msE) considered in this paper this criterion gives 



A > 0.24, 



(6) 



which is the same as for the Jupiter and a massless body, 
since mss/mj <C 1. We have found how good is this sim- 
ple analytical prediction in our case by performing direct 
orbit integration. Employing a simple N-body code, which 
uses Bulirsch-Stoer integrator and proved itself to be a con- 
venient way to obtain hi gh-accuracy solutions of the diffe r- 
ential equations (see e.g. Terquem and Papa loizou (2007)), 
we have calculated the planetary orbits in the region where 
the interior 3:2, 4:3 and 5:4 Jupiter mean-motion resonances 
are located. We have put the Super-Earth at the fixed loca- 
tion asE ~ 1, defining in this way our unit of distance, and 
the Jupiter initial semi-major axes have been varied in the 
range from 1.15 till 1.45 in these units. The initial eccentric- 
ities of both planets have been set to zero. Our results have 
been summarized in Table [T] where the number of conjunc- 
tions for each initial configuration that the system survives 
without suffering a close encounter is given in the column 3. 
When the Jupiter is located at the distance 1.24 or smaller 
we have observed the orbit crossing and scattering of the 
Super-Earth. Our result agrees ve ry well with the stability 
criterion given bv lGladmanl |1993). 

Hill stability criterion is useful for predicting close en- 
counters between planets but is not telling us much about 
other properties of motion a s for e xample whether or not the 
system is chaotic. IWisdoml (|l98Ch basing on the resonance 
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overlapping criterion found a boundary between chaotic and 
stable motions in the form 



A > 1.5ajm 



2/7 



(7) 



where m denotes mass ratio of the two biggest bodies in 
the system. The coefficient in front of the mass in the origi- 
nal paper was 2, but it has been imp roved later by numeri- 
cal simulations l|Duncan et al. |[l989l ) and we have adopted 
here the more recent value 1.5. The occurrence of the mean 
motion resonances is therefore very important for stability 
of the system. The width of the resonance regions for the 
circular r estrict ed three- body problem has b een derived by 
IWisdoml ljl980l l (see also iLecar et~aD ()200lh ). In Fig. [2] we 
have plotted the width of 3:2 and 4:3 i nterio r resonances for 
zero eccentricity given bv lLecar et al.l |2001) for the Jupiter 
(vertical dashed lines) . In this Figure the semi-major axis ra- 
tio is obtained by dividin g the semi-major axis of t he Jupiter 
by that of the small bodv. lWinter fc Murray! l|l997t ) have dis- 
cussed the various analytical models that have been used in 
the study of interior first order resonances. Analyzing sim- 
ple pendulum model they derived formula describing max- 
imum deviation of semi-major axis from the nominal value 
given by exact mean motion commensurability. Their re- 
sults for 3:2 and 4:3 resonances are presented also in Fig. [2] 
for comparison (solid lines). The theoretical prediction il- 
lustrated by the pendulum model (similar results were ob- 
tained using Hamiltonian approach) that for higher value of 
the eccentricity of the small body the first order resonance 
width is larger, has not been verified by direct integration 
ij Winter fc Murravl[l997h. The size of the resonant libration 
region derived by Winter fc Murray! (j 1997h performing nu- 
merical integrations of the full equ a tions of motion is quite 
similar to that given bv lLecar et al.l (|200ll ). The plot of max- 
imum libration widths has been supplemented with a few 
examples of the mean motion commensurability known in 
our Solar System. The dots denote Hilda and Thule groups 
of asteroids. This groups are particularly interesting because 
they are locked in the mean motion resonances with Jupiter. 
Note that asteroids are located exactly within regions pre- 
dicted by analysis of the restricted three-body problem. 

We have performed an analysis of the chaotic regions for 
our system (the Super-Earth and Jupiter around the Solar- 
type star) in the similar way it has been done for asteroid 
belt. First, let us evaluate the size of the chaotic region using 
equation (0. Substituting in this equation for m the mass 
ratio of Jupiter and Sun we have found that the small body 
(here the Super-Earth) should be outside the region of large 
scale chaos if 



A > 0.21aj. 



(8) 



According to this criterion the Super-Earth would be in the 
region of large scale chaos when the Jupiter's semi-major 
axis is smaller than 1.27. 

Next, we have examined the properties of the planetary 
orbits in the Hill stable region, which means for the Jupiter 
semi-major axis having value in the range between 1.24 and 
1.45. We have found two families of regular orbits separated 
from each other by the zone of the chaotic behaviour, which 
is located between 1.277 and 1.289. It should be pointed out 
here that the fact of the existence of the chaotic zone outside 
the region indicated by the resonance overlap criterion is not 
entirely unexpected. It should be kept in mind that criterion 
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Figure 2. Maximum libration regions in (semi-major axis ratio, 
eccentricity) space for the 3:2 (right) and 4:3 (left) resonances. 
The semi-major axis ratio has been calculated dividing the semi- 
major axis of Jupiter by the semi-major axis of the small body. 
The eccentricity is that of the small body. The dashed line is 
the resonance width for zero eccentricity and the solid line is the 
libration width derived using the pendulum approach. The dots 
correspond to the locations of asteroids in the 3:2 (Hilda group) 
and 4:3 (Thule group) mean motion resonances with Jupiter in 
the Solar System. We plot Thule asteroid and Hilda, Babylon, 
Normannia, Xanthomalitia, 2000VV39, 2001PR1 and 2002DL3 
from Hilda group using data from the Near-Earth Object Website 
hosted by The Jet Propulsion Laboratory. 



given by equation (0 is quite approximate and Wisdom's 
scaling law is strictly valid only in the asymptotic limit, 
when the integer p defini ng first order re sonance (p + 1) : p 
is much bigger than one (|Malhotr a||l996i). In our case of the 
Jupiter mass planet we are concerned with the first order 
resonances with small p. 

The families of the stable orbits can be nicely portrayed 
drawing the structure of the phase space trajectories in a 
plane on which the polar coordinates are (e, cf>) where e is 
the Super-Earth eccentricity and 4> ls the resonance angle 
defined as 



= 3As_b — 2Aj — zusi 



(9) 



Here Xse and Aj are the instantaneous mean longitudes 
of the Super-Earth and Jupiter respectively, and zuse is the 
longitudes of periastron of the Super-Earth. For given initial 
values of the eccentricity and resonance angle the motion of 
planets will follow the closed curve in this plane. In Fig. [3] 
we have shown three different phase space trajectories for 
the three different Jupiter locations at 1.27, 1.28 and 1.30. 
If the Jupiter is at 1.27, the Super-Earth has a regular very 
low eccentricity orbit (first panel), if it is placed at 1.30 the 
Super-Earth trajectory librates around the stable equilib- 
rium point in the 3:2 resonance zone (third panel), and in 
between (when the Jupiter is at 1.28) the planet experience 
a chaotic behaviour . Simi lar chaotic zone has been found by 
iHolman fc Murray! l|l99rj) in the asteroid belt. They have 
pointed out the importance of higher order resonances for 
the origin of chaos. In this context is worth mentioning the 
existence of 13:9 commensurability in the vicinity of this 
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Figure 3. The structure of phase space trajectories of the Super- 
Earth in three different planet configuration. From the top to the 
bottom panels, the Jupiter initial semi-major axis is 1.27, 1.28 and 
1.30 correspondingly. The distance from the origin is measured by 
the Super-Earth eccentricity, e and the polar angle is given by the 
resonance phase <j>. 



chaotic zone. This part of our investigations has been also 
summarized in Table [1] 

In the presence of gas giant, the dynamics of the Super- 
Earth with a mass of 5.5 M $ is similar to that of the massless 
body in the restricted three-body problem. This gives us 
an opportunity to determine some of the properties of its 
motion by means of ana lytical, well understood methods 
|Murrav fc Dermotufl999h . The extensive discussion of the 
classical three-body problem (Solar-type star, Jupiter and 
Super-Earth) brought us to the stage where we can predict 
stability and the character of the planetary orbits. At the 
same time we have prepared the well defined framework for 
our main study, which aim is to determine how this picture 
will be modified in a young planetary system in the presence 
of the gaseous protoplanetary disc in which the planets are 
embedded. 



4 THREE-BODY PROBLEM WITH A 
GASEOUS DISC 

4.1 Description of the performed simulations 

The system described in the previous Section has been con- 
sidered again, but this time the planets are embedded in 
a gaseous disc. We investigate the orbital evolution of the 
Super-Earth and Jupiter in the same region as before tak- 
ing into account the disc-planet interaction. The disc has a 
constant aspect ratio — = 0.05. It extends from 0.33 till 
4 (the unit of distance is the initial semi-major axis of the 
Super-Earth as it has been chosen in the previous Section). 
The kinematic viscosity is constant in space and its value 
is v — 10~ 5 in units of a 2 SE (GM, /a% E ) 1 ^ 2 , where asE is 
the initial semi-major axis of the Super-Earth and G is the 
gravitational constant. This value, with a disc aspect ratio 
of 0.05 should correspond to a = 4- 10~ 3 at the distance 1 in 
our units in the standard a-parametrization for the viscosity 
in the disc. The planets has been placed in the fiat part of 
the surface density profile with E = 2 • 10 3 (5.2au/as£;) 2 kg 
m~ 2 : the standard value attributed to the minimum mass 
solar nebula at 5.2 au when clse = 5.2 au. The initial ec- 
centricity of both planets is set to zero. We perform nu- 
merical simulations using hydrodynamical code NIRVANA, 
for details of the n umerical scheme and code adopted, see 
iNelson et ail (|200Ch . We have used the resolution 384 x 512 
grid-cells in the radial and the azimuthal direction respec- 
tively. We are aware of the fact that adopted resolution 
might be not high enough to simulate properly the evolution 
of the Super-Earth in general, because the corotation region 
around planet is covered only by a few grid cells. However, 
we think that adopted resolution is a good compromise be- 
tween the accuracy of the simulations and reasonable com- 
putational time. The computational domain is a complete 
ring where the azimuthal coordinate ip — corresponds to 
<p — 2n. We put the open boundary conditions in the ra- 
dial direction. The potential is softened with the parameter 
b — 0.04r in such a way that the 1/r law has been changed 
into l/Vr 2 + b 2 . 

As in the case described in previous Section, where the 
disc was absent we have investigated 5.5M0 and lMj plan- 
ets with the initial positions specified in Table[T] The Super- 
Earth is located initially at the distance 1 in each simulation. 
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Table 1. The initial separation of planets in Super-Earth units (in Jupiter units) and their final configurations. The double 
line denotes the Hill stability boundary. 



Initial Super-Earth Initial Jupiter Final outcome of the evolution Final outcome of the evolution 

location location without disc with disc 



1.0 


(0.870) 


1.15 


(1.0) 


scattering (before first conjunction) 


scattering 


1.0 


(0.862) 


1.16 


(1.0) 


scattering (before first conjunction) 


scattering 


1.0 


(0.855) 


1.17 


(1.0) 


scattering (before first conjunction) 


scattering 


1.0 


(0.847) 


1.18 


(1.0) 


scattering (after 39 conjunctions) 


scattering 


1.0 


(0.840) 


1.19 


(1.0) 


scattering (after 20 conjunctions) 


scattering 


1.0 


(0.833) 


1.20 


(1.0) 


scattering (after 2 conjunctions) 


4:3 


1.0 


(0.826) 


1.21 


(1.0) 


scattering (after 15 conjunctions) 


scattering 


1.0 


(0.820) 


1.22 


(1.0) 


scattering (after 9 conjunctions) 


4:3 


1.0 


(0.813) 


1.23 


(1.0) 


scattering (after 107 conjunctions) 


4:3 


1.0 


(0.806) 


1.24 


(1.0) 


scattering (after 5 conjunctions) 


4:3 


1.0 


(0.800) 


1.25 


(1.0) 


survival (more than 100,000 conjunctions) 


4:3 


1.0 


(0.794) 


1.26 


(1.0) 


survival (more than 100,000 conjunctions) 


3:2, scattering 


1.0 


(0.787) 


1.27 


(1.0) 


survival (more than 100,000 conjunctions) 


3:2, scattering 


1.0 


(0.781) 


1.28 


(1.0) 


survival (more than 100,000 conjunctions) 
chaotic orbit 


3:2 


1.0 


(0.769) 


1.30 


(1.0) 


survival (more than 100,000 conjunctions) 


3:2 


1.0 


(0.741) 


1.35 


(1.0) 


survival (more than 100,000 conjunctions) 


3:2 


1.0 


(0.690) 


1.45 


(1.0) 


survival (more than 100,000 conjunctions) 


3:2 



The Jupiter mass planet is located on the external orbit at 
the different locations. Such configuration allows to attain 
convergent migration (Fig. [TJ. 

4.2 Super-Earth in mean motion resonance with 
Jupiter 

In agreement with the theoretical predictions (Fig.[TJ for our 
particular choice of the system parameters, the Jovian mass 
planet migrates faster than the Super-Earth, and thus plan- 
ets approach each other and may become eventually cap- 
tured in a mean motion resonance. In order to investigate 
how an outcome of such evolution depends on the initial 
configurations of planets we have performed our simulations 
with a wide range of the planetary separations (see Table [IJ. 
We are particularly interested in the conditions for which the 
first order mean motion resonances p + 1 : p for p ^ 2 are at- 
tained (here p is an integer). The 2:1 commensurability was 
not taken into account because of the computational time 
constraints. The resonances located closer to the Jupiter like 
5:4, 6:5 and with higher p are not possible, because for such 
small separations the system become unstable and Super- 
Earth is scattered from the disc. We have found that two 
outcomes are possible, namely the Super-Earth is ejected 
from the disc or planets become locked in 3:2 or 4:3 mean 
motion resonances. 

Our results has been summarized in Table [1] column 4 
and on the right hand side of the Fig. [4] where the most 



likely outcome of the planet evolution is given for every ini- 
tial semi-major axis ratios. In Fig. [4] we show also two exam- 
ples of the resonant trapping occurred for the relative planet 
separation 0.45 (upper curve) and 0.24 (lower curve). In the 
first case, the differential migration brought planets into 3:2 
commensurability and in the second one into 4:3. The semi- 
major axis ratio of planets librates exactly within regions of 
the resonance width marked by the solid, horizontal lines in 
Fig|4] We plot there the resonance widths for (p+1) : p com- 
mensurabilities with 2 p ^ 4 in a case of circular orbits 
in the restricted three b ody problem where the secondar y 
is a Jovian mass planet l| Wisdom! 1 1980l : iLecar et al.ll200lh . 
The width of 3:2 and 4:3 resonances have already been illus- 
trated in Fig. [2] and discussed in the previous Section. As we 
can see the 4 : 3 and 5 : 4 mean motion resonances partially 
overlap each other, which might have a significant effect on 
the stability of the system. Indeed planets located initially in 
this region were immediately scattered from the disc. Both 
commensurabilities shown in Fig. 3] are stable for the whole 
time of our simulations. We have followed their evolution 
for not longer than 2500 orbits (the unit of time is defined 
as the orbital period on the initial orbit of the Super-Earth 
divided by 2n) because after that time the surface density 
is significantly reduced due to viscous evolution of the disc. 
However, not all planetary pairs initially locked in the mean 
motion resonance survived further evolution (always limited 
to 2500 orbits). In the case of 3:2 commensurability the li- 
bration width increases with time. The amplitude of the 
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scattering 
4;3 resonance 



> scattering 



1 000 1 500 

time (orbits) 



2500 



Figure 4. The evolution of the ratio of semi-major axes for the 
Jupiter and Super-Earth embedded in a gaseous disc. In the case 
of upper curve the initial semi-major axis ratio is 1.45 and planets 
became locked into 3:2 resonance. For the lower curve, the initial 
semi-major axis ratio is 1.24 and the attained commensurability is 4:3. 
The solid horizontal lines are the 3:2, 4:3 and 5:4 resonance widths 
in the circular three-body problem. Note that 4:3 and 5:4 resonances 
partially overlap each other. 



oscillations is bigger for planets with a separation smaller 
than the lower boundary of the 3:2 resonance region. So in 
fact in the case of Jupiter placed at 1.26 and 1.27 the am- 
plitude grows and eventually the Super-Earth is scattered 
from the disc. We have illustrated it for the Jupiter at 1.26 
in Fig. [5] The direct reason of this scattering is the excita- 
tion of the Super-Earth eccentricity to the high value, which 
favours close encounters and ends with the ejection of the 
planet. 

We have also investigated eccentricity evolution of the 
Super-Earth in few other cases. Fig. [6] shows eccentricity of 
the smaller planet versus semi-major axis ratio. The dashed, 
vertical lines denote the resonances width for circular case, 
the solid lines instead the maxi mum libration wid t h ob - 
tained from the formula given by IWinter fc Murray! (I1997T I 
in the pendulum approach in the case of the classical re- 
stricted three body problem. We have already introduced 
these quantities in the previous Section. The eccentricity of 
the Super-Earth locked in 3 : 2 commensurability increases 
and at the end of simulation reaches value about 0.5. Eccen- 
tricity of planets approaching 3 : 2 mean motion resonance 
from higher value of semi-major axis exceed value of 0.5. 
For the planets attaining commensurability from the lower 
value of the semi-major axis ratio than needed for the exact 
3:2 resonance, the oscillations around exact position of the 
resonance increase and eventually Super-Earth is scattered 
from the disc (Fig. [5]). 

The Super-Earth locked in 4:3 mean motion resonance 
in all simulations approach eccentricity value around 0.3. 
This is a lower value of eccentricities than in the case of 
3:2 commensurability. The amplitude libration around the 
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1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 
semi-major axis ratio 

Figure 5. The Super-Earth eccentricity excitation for the case of 
Jupiter located at 1.26. The amplitude of the semi-major axis ra- 
tio oscillations around the position of the 3:2 resonance increases 
and eventually the Super-Earth is scattered from the disc. The 
dashed line is the zero eccentricity resonance width and the solid 
line is the libration width obtained in the pendulum approxima- 
tion for 3:2 resonance (the same as in Fig [2j . 

position of the exact resonance is smaller than for 3:2 and 
there is no trend in the amplitude to increase. 

Till now we have concentrated on the orbital elements 
evolution occurring due to the tidal planet disc interaction. 
Another effect of this interaction in the case of the giant 
planet is the gap opening process, which changes the sur- 
face density profile and may affect the Super-Earth migra- 



8 



E. Podlewska and E. Szuszkiewicz 



o 




1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 



semi-major axis ratio 

Figure 6. Four examples of the Super Earth eccentricity evo- 
lution. For the 4:3 commensurability (left) we are presenting the 
eccentricity as a function of semi-major axis ratio for the plan- 
ets with their initial separation 0.24 (gray colour) and 0.22 (black 
colour) and for the 3:2 commensurability (right) the initial planet 
separation reads 0.45 (gray colour) and 0.27 (black colour). The 
resonance widths are the same as in Fig. [2] 



0.0008 



0.0007 - 



0.0006 - 



a) 

"2 0.0004 - 
o 

1 0.0003 - 
to 

0.0002 - 



0.0001 - 








0.5 1 1.5 2 2.5 3 3.5 4 
radius 

Figure 8. The disc profile and planet locations at the time of 4:3 
resonant capture. The initial position of Jupiter is 1.24. 



and ongoing process of the gap opening in the disc leads to 
the resonant configurations for a range of the initial planet 
locations. 
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Figure 7. The disc profile and planet locations at the time of 3:2 
resonant capture. The initial position of Jupiter is 1.30. 



tion. In Fig. [7] we show the surface density profile together 
with the positions of planets at the moment when they ap- 
proach 3 : 2 mean motion resonance. In this simulation the 
Jupiter is placed at 1.30 which is close to the location of 
the nominal resonant semi-major axes. After about 628 or- 
bits the resonant behaviour is clearly seen in both resonance 
angles. The gap is still not entirely opened. The Jupiter is 
located in it and the Super-Earth is right at the edge of the 
gap. In Fig. [8] we plot the same as in Fig. [7] but for 4 : 3 
commensurability. Here the Jupiter started at 1.24 and the 
resonant capture happened after about 1178 orbits. It can 
be seen that at the time of locking in 4 : 3 mean motion 
resonance the Super-Earth orbits in lower surface density 
region in the steep gradient of surface density at the wall of 
the gap. 

The interplay between the migration of both planets 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the evolution of Super- 
Earth and Jupiter mass planets embedded in gaseous pro- 
toplanetary discs where terrestrial planet is moving on the 
internal orbit and gas giant on the external one. We have 
found that according to the simple theoretical prediction dif- 
ferent rates of migration lead to the planet capture in the 
mean motion resonances. For the system considered here the 
only possible first order commensurabilities are with p ^ 3. 
For higher p the resonance locations will be in the Hill un- 
stable zone and in the regions where the resonances begin 
to overlap therefore the system will be unstable. 

The essential feature of our scenario or differently, the 
reason for the occurrence of the mean motion resonances in 
our system, is the convergent migration of planets embed- 
ded in the protoplanetary discs. We have mentioned already 
in Section [2] that for the given masses and disc profile we 
can control the convergence of the migration by the vis- 
cosity in the parent disc, at least for the standard range 
of parameters, as in the case considered here. In general, 
the situation may be more complex. First of all, it has 
been found that the torque acting on the low-mass plan- 
ets might be considerably smaller than the analytic lin- 
ear estimate and its va lue does depend on the viscosity 
(see iMasset et alj {2006) for the discussion). Moreover, it 
has been also proposed that even the planets as small as 
Super-Earths can open a gap in the dis c and slow down 
their orbital migration l|Matsumura et all ([2007) and refer- 
ences therein) . Finally, combining the equation @ with the 
gap openi ng formula suggested fo r a disc with a very small 
viscosity \hin fc Papaloizoul Il993r i . ICrida et"all (|2006t ) con- 
structed the generalized criterion which is more restrictive 
than equation Q. It is than required to consider in detail 
the disc profile, masses of the planets and their locations 
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in the disc in order to determine whether or not we might 
expect the convergent migration. 

Our results for the case of the mature planetary systems 
when the disc is absent are consistent with those predicted 
analytically for the classical planetary three-body problem. 
From the direct orbit integration we have found the Hill sta- 
bility region which occurs when the Jupiter and Super-Earth 
orbit relative separation is 0.24 or bigger, which is exactly 
what we get from the formula (B). The chaotic orbits have 
been identified in the configurations with the initial semi- 
major axis ratios in between 1.277 and 1.289. This chaotic 
zone might be connected with 13:9 resonance, which par- 
tially overlap with the 3:2 commensurability. 

As has been summarized in Table Q] the planets em- 
bedded in the disc show the same behaviour as planets in 
matured planetary system (i.e. where the disc is absent) for 
the relative separation between planets less than 0.20. So 
for these separations the presence of the gas in the system 
does not change the outcome of the evolutio n. The planets 
are unstable in the Hill sense (lGladmanlll993l ) and also they 
are in the region with the resonance overlapping dWisdoml 
1980). The first difference between young and mature plane- 
tary systems has been found for the relative separation 0.20. 
The Super-Earth has been trapped by Jupiter in 4:3 mean- 
motion resonance very quickly preventing in this way the 
scattering. This was not the case for 0.21 and the Super- 
Earth has been ejected from the disc. For the even higher 
separations from 0.22 till 0.24 which are also Hill unstable 
we have obtained the 4:3 resonance. The last 4:3 commen- 
surability we have got for the Super-Earth separated from 
Jupiter by 0.25. Next two configurations with se parations 
0.26 a nd 0.27, located inside the chaotic zone of IWisdorr] 
(1980) passed through unstable 3:2 resonance. All other in- 
vestigated by us structures, namely for the separations 0.28- 
0.45 have been locked in 3:2 commensurability. 

The presence of the disc and its interactions with the 
planets might be a key to the attaining resonant configu- 
rations like those observed between Jupiter and groups of 
asteroids in our system. The fate of the Super Earth in the 
case of the presence of Jupiter mass planet on the external 
orbit is determined by separation of the planets and their 
eccentricity evolution. The resonant structures with 3:2 and 
4:3 commensurabilities are easy to attain but their stability 
requires further studies. Our scenario indicates the possi- 
bility that terrestrial planets and giant planets coexist in 
relatively close proximity. It has an interesting implication 
for the astrobiological studies. For example, if we consider 
a system, which contains a giant planet in or near its hab- 
itable zone, taking into account a simple stability analysis 
as in Section [3] we might conclude that there is no dynam- 
ical room for terrestrial planets in the habitable zone, close 
to the Jupiter mass planet. However, as it has been shown 
here, if we allow for the possibility of mean motion resonance 
locks among planets, this is no longer necessarily the case. 
This hypothesis can be verified by the discoveries of Jovian 
type planets in the habitable zones, similar to the low mass 
gas pl anet which has been found recently by iFischer et al.l 
(2007) in 55 Cancri system, so the intensive search for the 
low mass planets which may be present close to the gas gi- 
ants is more than justified. The discovery of the resonant 
configurations would provide constraints also on the migra- 
tion processes themselves. If the future observations will find 



primarily smaller planets locked in the interior mean motion 
resonances and not those locked in the exterior ones as in 
iThommesI (I2005D . or other way round, than we will get im- 
portant information on the characteristics of the protoplan- 
etary disc in which these planets form. Finally, the occur- 
rence of mean motion resonances may increase the chance 
of a detection of the terrestrial planets. The best example 
is provided by the transit timing measurements, which may 
detect additional planets (which are as small as the Super- 
Earth considered here) in the system via their gravitational 
interaction with the transiting Jupiter mass planet, taking 
advantage of the resonance induced by migration. 
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